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We present a theoretical study of ghost imaging based on correlated beams arising from parametric 
down-conversion, and which uses balanced homodyne detection to measure both the signal and idler 
fields. We analytically show that the signal-idler correlations contain the full amplitude and phase 
information about an object located in the signal path, both in the near-field and the far-field case. 
To this end we discuss how to optimize the optical setups in the two imaging paths, including the 
crucial point regarding how to engineer the phase of the idler local oscillator as to observe the desired 
orthogonal quadrature components of the image. We point out an inherent link between the far-field 
bandwidth and the near-field resolution of the reproduced image, determined by the bandwidth of 
the source of the correlated beams. However, we show how to circumvent this limitation by using 
a spatial averaging technique which dramatically improves the imaging bandwidth of the far-field 
correlations as well as speeds up the convergence rate. The results are backed up by numerical 
simulations taking into account the finite size and duration of the pump pulse. 

PACS numbers: 42.50.Dv,42.50.-p,42.50.Ar,42.30.Wb 



I. INTRODUCTION 

Ghost imaging relies on the spatial correlation between 
two beams created W, e^g.-parametric down conversion 

(pdc) [aaiaiEBiaiiaiiiiiiEQia Each 

of the correlated beams are sent through a distinct imag- 
ing system called the test arm and the reference arm. In 
the test arm an object is placed and the image of the ob- 
ject is then recreated from the spatial correlation function 
between the test and reference arm. A basic requisite of 
the ghost-imaging schemes is that by solely adjusting the 
reference arm setup and varying the reference point-like 
detector position it should be possible to retrieve spatial 
information about the object, such as the object image 
(near field) and the object diffraction pattern (far field). 

Traditionally the studies devoted to ghost imaging 
have considered PDC in the low-gain regime, where the 
conversion rate of the pump photons to a pair of en- 
tangled signal-idler photons is low enough for the detec- 
tor to resolve them one at a time. This leads to the 
so-called two-photon imaging schemes first investigated 
by Klyshko [J, The recreation of the object is then 
based on coincidence counts between the test and refer- 
ence arm, and the working principles have been demon- 
strated experimentally [3, U, |5j . Our group has focused 
on generalizing the governing PDC theory to the macro- 
scopic, high-gain regime where the number of photons per 
mode is large pj| [Tjl HiA |2(j ■ Moreover, we general- 
ized the theory behind the two-photon i mag in g schemes 
H II B II II 13 to the high-gain regime El |H 13 and 
also investigated the quantum properties of the signal- 
idler correlations in that case [Til lisl H^. |20| . While co- 
incidence detection is used for low gain, in the high-gain 
regime the object information is extracted from measur- 
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ing pixel-by-pixel the signal-idler intensities and from this 
forming the correlations. 

Here we consider a setup where each of the PDC beams 
are measured with balanced homodyne detection by over- 
lapping them with local oscillator (LO) fields. A simi- 
lar scheme was studied previously by our group in the 
context of entanglement of the signal-idler beams from 
PDC where it was used to show analytically that 
the entanglement is complete since it encompasses both 
the amplitudes and the phases of signal and idler. The 
initial motivation for using a homodyne scheme for quan- 
tum imaging came from the need to circumvent the prob- 
lems related to information visibility in the macroscopic 
regime. Specifically, when intensity detection is per- 
formed a homogeneous background term is present in the 
measured correlation function [H 13 H3- This term, 
which can be rather large, does not contain any infor- 
mation about the object and lowers the image visibil- 
ity. Instead, by using homodyne detection the signal- 
idler correlation becomes second order instead of fourth 
order, and hence this background term is absent. An- 
other advantage of homodyne detection is that arbitrary 
quadrature components of the test and reference beams 
can be measured, which means that the homodyne de- 
tection scheme allows for both amplitude and phase mea- 
surements of the object. We will show detailed analytical 
calculations which demonstrate that this is indeed pos- 
sible, both for the object image (near field) as well as 
its diffraction pattern (far field). It is possible to recon- 
struct even a pure phase object when a bucket detector is 
used in the test arm (in contrast to when intensity mea- 
surements are done, as it was shown for the coincidence 
counting case 01)- 

We will present a technique that implements an aver- 
age over the test detector position, and this turns out 
to strongly improve the imaging bandwidth of the sys- 
tem with respect to when a fixed test detector position is 
used, where the imaging bandwidth is limited to the far- 
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field bandwidth of PDC. Additionally the method sub- 
stantially speeds up the convergence rate of the far-field 
correlations. This implies that even complex diffraction 
patterns with high-frequency Fourier components can be 
reconstructed in a low number of pump pulses. 

The suggested homodyne detection scheme is fairly 
complicated to implement experimentally, since it in- 
volves carrying out two independent homodyne measure- 
ments of fields that are spatially and temporally multi- 
mode (with pulse durations on the order of ps). A dual 
homodyne measurement has been done in the case of spa- 
tially single- mode light pi I22T] . but even a single homo- 
dyne measurement of spatially multimode light has not 
previously been investigated in details to our knowledge, 
not even in the continuous wave regime where the field is 
temporally single-mode. Homodyne experiments of spa- 
tially single-mode but temporally multimode light have 
pointed out the importance of a proper overlap between 
the LO and the field |23J, |24|, and alternative ways of 
producing the LO has been implemented [25J leading to 
a better overlap. An improper overlap leads to a degrad- 
ing in, e.g., the squeezing as it effectively corresponds 
to a drop in the quantum efficiency [2(|. In our case 
these problems are presumably less severe since the quan- 
tum efficiency does not play an important role. In fact, 
we will present numerical simulations that confirm that 
homodyne measurement of spatio-temporal multimode 
light can be done successfully, with the strongest techni- 
cal restriction concerning the engineering of the phases of 
the local oscillators. It is worth to mention that the ho- 
modyne measurement protocol (no background term in 
the correlations) in combination with the spatial averag- 
ing technique (increased spatial bandwidth) could open 
for the possibility of using an OPO below threshold for 
imaging. This would simplify a possible experimental 
implementation. 

In Sec. [n] the model for PDC is presented, and the 
the optical setup of the imaging system is described in 
Sec. IIIII In Sec. IIVI we discuss the case where point-like 
detectors are used in the test arm, while in Sec. [3 the 
bucket detector case is discussed, and show analytical 
results in the stationary and plane-wave pump approxi- 
mation concerning the retrieval of the image information 
of both the far and the near field, as well as numerical 
results that include the Gaussian profile of the pump. 
Besides confirming the analytical results these serve as 
examples for discussion and demonstration of the imag- 
ing performances of the system. In Sec. I VII we draw the 
conclusions. The appendices discuss a quadratic expan- 
sion of the PDC gain phase (App. [SJ, and how to treat 
the temporal part of the analytical correlations (App.[5J|. 



II. THE MODEL 

The starting point is a general model describing the 
three-wave quantum interaction of PDC inside the non- 
linear crystal |l9|, |20j , which includes the effects of finite 



size and duration of the pump and the effects of spatio- 
temporal walk-off and group-velocity dispersion. 

We consider a uniaxial nonlinear crystal cut for 
type II phase matching. The injected beam at the fre- 
quency loq (pump beam) is sent into the crystal in the z 
direction. Inside the crystal the pump photons may then 
down-convert to sets of photons at the frequencies u>\ 
(signal photon, ordinarily polarized) and L02 (idler pho- 
ton, extraordinarily polarized). The classical equations 
governing PDC are in the paraxial and slowly varying 
envelope approximation 

CjAj = <iA Ate^ A " z , j,l = l,2, 3+I (la) 
C A = -aA 1 A 2 e iAoZ . (lb) 

Cj are operators describing linear propagation in the 
medium, including walk-off as well as 2nd order disper- 
sion effects in the spatial and temporal domains 

8 ,,d ik" d 2 d i „ 2 

^ = ^ + ^ + ^-^-2^- V - < 2 > 

kj = rijUjj/c is the wave number of field j inside the 
crystal, and the primes denote derivatives with respect 
to u taken at the carrier frequency tOj : the two terms de- 
scribe temporal walk-off (kj = dkj/dcj\ u=ulj ) and group 
velocity dispersion (k'J = d 2 kj /du?\ u =Wj)- Dealing for 
simplicity with uniaxial crystals, we have assumed that 
the walk-off direction of the extraordinary waves is along 
the re-axis. The walk-off angle pj = dk/dq x \^ = o deter- 
mines the energy propagation direction of wave j with 
respect to the pump axis z. The effect of diffraction is 
described in the framework of the paraxial approxima- 
tion through the Laplacian = d 2 /dx 2 + d 2 jdy 2 with 
x = (x, y) spanning the transverse plane. 

The right hand sides of Eq. Q give the nonlinear in- 
teraction in the crystal, the strength of the nonlinearity 

being governed by a that is proportional to the effec- 

(2) 

tive second order susceptibility \ c g ■ ^0 = k\ + k^ — ko 
denotes the collinear phase-mismatch between the three 
waves along the z-axis. 

These classical equations can be converted into a set of 
operator equations describing the evolution of the quan- 
tum mechanical operators. This conversion is formally 
done by substituting Aj(x, z, t) — ► dj(x,z,t), j = 1,2. 
a,j(x, z, t) are boson operators obeying the following com- 
mutator relations at a given z 

[ai(z,x,t),a!j(z,x?,t')] = 5 rj 8(x - x)8(t - t'), (3) 

j = 1,2, while all other combinations commute. With 
this choice the expectation value of the intensity aJa gives 
the number of photons per time per area. The pump Aq 
is treated classically, and furthermore we adopt the para- 
metric approximation where the pump is taken as unde- 
pleted. It is more convenient for the subsequent analysis 
to present the equations in Fourier space according to the 
following transformation 

aj(z,q,n) =J§J ^ aj (z,x,t)e-^+^, (4) 
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j = 1,2. The signal equation i|la|) is then 
C 1 (z,q,n)a 1 (z,q,n) = ae- iA °* f ^ 



xAq(z, q — if, SI — Sl')a 2 (z, ■ 



-q , 



dsi 1 
-si'), 



(5) 



while the idler equation can be found by exchanging sub- 
scripts 1 ^> 2. The Fourier version of the linear propa- 
gation operator is defined as Cj(z, q, SI) = — i5j(q, SI), 
where the advantage of being in the Fourier space is that 
the derivatives with respect to t and x become constants 



1 



s 3 (q,n) = k l J n + -k';n 2 + Pj q, 



i 

2k j 



(6) 



The classical pump is taken as being Gaussian in both 
time and space and can be expressed as 

A (z = 0,x,t) = {2irf' 2 A p e-W 2 '<- t2/T \ (7a) 
2V2Ap -|g|'/ggg-o a /Jh>g 



A (z = 0,g,O) 



(7b) 



in the real and Fourier space, respectively. Here Wq and 
To are the pump waist and duration time, while A p is the 
pump amplitude. The bandwidths in Fourier space are 
given by Sqo = 2/wq and Suj® = 2/tq. 

In the stationary and plane- wave pump approximation 
(SPWPA) the pump is taken as translationally invariant 
as well as continuous wave. Thus, wq and To tend to 
infinity so we have A (z,q,Q) -» (2nf/ 2 A p 5(q)5(Sl). 
Under this condition Eqs. (|3J) can be solved analyti- 
cally. The unitary input-output transformations relat- 
ing the field operators at the output face of the crystal 
a° ut (g, St) = a,j(z = l c ,q,Sl) to those at the input face 
a™(q, St) = aj(z = 0, q, St) then take the following form 

a?* ft SI) = UMn)a?(q,n) 

+ Vi(q, Siy^i-q-, -SI), i^j = 1, 2. (8) 
The signal gain functions are 



U 1 {q,Sl)= e *Di2(q,n)ij2 cosh(T 12 (q,Sl)l c ) 



A 12 (g,n) 

2r 12 (g,fi) 



sinh(r 12 (g,fi)/ c ) 



Vi(q,Sl) 



r iDi ■>(Q,U)lrl2 sinh(r 12 (<f,Sl)i c ) 

e r 12 (q.n) 



(9a) 
(9b) 



For the idler similar gain functions are found by ex- 
changing indices 1^2. We have in Eq. introduced 



a v = crAp as well as 



Dij(q,Sl) = 5 i (q,Q)-5 j (-q,-Sl.)-A c 



^-[A„(9,0)/2P, 



(10a) 
(10b) 
(10c) 



Ty {q, SI) 

Ay(£ft) = A + 8 t (q,Sl) + S :j (-q,-Sl). 

It is important to note that the gain functions satisfy the 
following unitarity conditions 

\U 3 (q^)\ 2 -\V {m\ 2 = 1, (i = 1, 2) (11a) 
Ut(q, Sl)V 2 (-q, -SI) = U 2 (-q, -Sl)Vx(q, SI), (lib) 



which guarantee the conservation of the free-field com- 
mutation relations © after propagation. 

The bandwidths of emission of the gain functions © 
in the spatial and temporal frequency domain are 



ql = 2[Z c (l/fc 1 + l/fc 2 )]-\ 
Sl = - fc^c) -1 , 



(12a) 
(12b) 



where the typical variation scales of the signal-idler fields 
in space and time are 



Xcoh = l/<70, T co h = 1/SIq- 



(13) 



They can be identified with the coherence length and 
coherence time, respectively, of the fields. 



III. THE SYSTEM SETUP 

Ghost imaging is characterized by its two-arm config- 
uration, with an unknown object placed in the test arm. 
The information about the object is retrieved from the 
cross-correlations of the fields recorded in the test and 
reference arms as a function of the reference arm pixel 
position, see Fig.^ By simply changing the optical setup 
in the reference arm, information about both the object 
image (near field) and the object diffraction pattern (far 
field, i.e. the Fourier transform) can be obtained. There 
are two main motivations for using such an imaging con- 
figuration. Firstly, this configuration makes it possible 
to do coherent imaging even if each of the two fields are 
spatially incoherent (and thus recording the far-field spa- 
tial distribution in the test arm does not give any infor- 
mation about the object diffraction pattern). Secondly, 
the configuration allows for a simple detection protocol in 
the test arm even without measuring spatial information, 
and still the spatial information can be retrieved from the 
correlations. Thus a simple bucket detector setup can be 
used that collects all photons, or alternatively using a 
single pixel detector (point-like detector). This is advan- 
tageous when the object is located in a environment that 
is hard to access making it difficult to place an array of 
detectors after the object, or if the stability of the test 
arm is an issue making a simple setup crucial. 

In this publication we will focus mainly on the point- 
like detector in the signal arm for two reasons. First 
of all, in the high-gain regime there are many photons 
per mode implying that a single-pixel measurement carry 
enough information to perform the correlation. In con- 
trast, in the low-gain regime the probability of having a 
photon after the object is so low that a bucket detector 
is sometimes favored in order to get coincidences within 
a reasonable time. Secondly, having spatial resolution in 
the signal arm allows to perform averages over the pix- 
els of the test arm detector (in this case it is necessary 
to have an array of pixels in the test arm), which as we 
will show vastly speeds up the convergence rate as well 
as improves the imaging bandwidth. For completeness 
we also treat the bucket detector case in Sec. El 
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FIG. 1: The system setup. Sketch (a) shows the general setup 
where a point-like detector (or a pixel array) is used in the test 
(signal) arm. Sketch (b) shows the case of a telescope setup 
in the reference (idler) arm. Sketch (c) shows the test arm 
setup when bucket detectors are used for the measurements 
there. See the text for more details. L: lens of focal length /, 
PBS: polarized beam splitter, BS: beam splitter, P: crystal 
plane, P{: signal measurement plane, P^\ idler measurement 
plane. 



We consider the case where the signal and idler input 
fields are in the vacuum state. Thus, in the SPWPA we 
have the following correlations at the crystal exit [Ifl 12(1 ] 



= 5 jk 5(q- cf)5(n - fi')| W> ^)| 2 , 3 = 1,2 (14a) 

= 5{q + ?)5(Sl + n')U x {q, Q)V 2 (-q, -fi). (14b) 

The type II phase matching conditions of the crystal 
ensure that the signal and idler have orthogonal polariza- 
tions and therefore at the exit of the crystal they can be 
separated by a polarizing beam splitter (PBS). We will 
neglect the distance between the crystal and the PBS. 
The fields Cj (x, t) at the measurement planes Pj are con- 
nected to the output fields at the crystal exit plane P by 
a Fresnel transformation formally written as 

Cj(xj,t) = dx*h j (x j ,x*)a l j nt {x',t) + V j (x j ), (15) 



where hj(xj,x') are kernels related to the optical path 
from P to P'y Vj(xj) accounts for losses in the imag- 
ing system that are linearly proportional to vacuum field 
operators and therefore uncorrelated to the PDC fields. 

The object is described by a transmission function 
TobjO?), and is set in the signal arm. Besides the object 
the signal arm is set in the so-called f-f imaging scheme, 



consisting of a lens L with focal length / located at a 
distance / from both the crystal exit and from the detec- 
tion plane P[. If we apply the unfolded Klyshko picture 
it is obvious that the required position of the object 
changes depending on the type of detector used in the 
test arm. Namely, when a point-like detector is used the 
detector acts like a point-like source which through the 
f-f lens system becomes plane wave source at the crys- 
tal, and therefore the optimal object placement is at P, 
as shown in Fig. QJi. Conversely, if a bucket detector is 
chosen it acts like a plane wave source at P{ and through 
the f-f lens system it changes to a point-like source at P. 
Thus, if in this case we put the object at P we will not 
get any information about it, while if we put the object 
at P[ it is instead effectively illuminated by the plane 
wave source. Thus, in the bucket case the object must 
be placed at P{, as shown in Fig. Consequently, the 
signal Fresnel kernels are 

htfutf) = (iXxfy'e-^^T^af), (16a) 
h^xux") = (zAr/rV^^VobjOrx), (16b) 

for the point-like and the bucket detector case, respec- 
tively, and the superscript "b" denotes bucket. Ai is the 
wavelength of the signal field. Note that the lens in the 
signal arm could actually be removed and substituted by 
a propagation over a distance long enough to be in the 
Frauenhofer regime. However, the phase of the signal 
field could be difficult to control in that case. 

We will fix the setup of the signal arm once and for all. 
We may then adjust the idler arm to retrieve the desired 
kind of information through the signal-idler correlations. 
The idler lens system also consists of a lens L of focal 
length /, as shown in Fig.^i. In one setup d± and d 2 are 
identical and both equal to the lens focal length, d\ = 
c?2 = /, i.e. an f-f setup as also used in the signal arm. 
With this setup the correlations in the point-like detector 
case give information about the object far field, while 
in the bucket detector case they give information about 
the object near field. In the other setup, instead, these 
distances obey the thin lens law l/di + l/d-z = 1/f- This 
means that the object near field can be retrieved from 
the correlations in the point-like detector case, while the 
object far field can be retrieved in the bucket detector 
case. For simplicity we may keep di = d% = 2f and we 
denote this the 2f-2f setup. We will also consider the so- 
called telescope case, see Fig.^D, where the 2f-2f setup is 
repeated twice adding an extra lens, but the system still 
obeys the thin lens law: therefore the object information 
can be retrieved as in the 2f-2f case. For these three 
chosen setups in the idler plane, the idler field at P 2 is 
then given by the Fresnel transformation (|15fl with the 
kernels 



h 2 j(x,x") 

h2,2f(x,x') 
h2,T{x,x') 



(lAa/)-^-*^*-* , 
5(x — x), 



(17a) 

(17b) 
(17c) 
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for the f-f setup, the 2f-2f setup, and the telescope setup, 
respectively. We choose to work with the telescope setup 
with respect to the 2f-2f setup because of the direct trans- 
formation from the image plane P to the measurement 
plane P 2 [due to the 5(x' — x) type correlation], as well 
as the lack of the phase factor. 

The fields at P'j are measured using balanced homo- 
dyne detection schemes, so Cj(xj, t) is mixed with an LO 
aj(xj,t) on a 50/50 beam splitter, and after the beam 
splitter the two fields are 

Cj (x, t) ± Oj (x, t) 

Cj,±{X,t)= -y= . (18) 

The LO is treated as a classical coherent field, therefore 
by measuring the field intensities at the photo-detectors 
and subtracting them we obtain 

Zj(x,t) = c' j + (x,t)cj i+ (S,t) - & jt _(S,t)cj-(S,t) 

— cj (x, t)a* (x, t) + Cj (x, t)ctj (x, t) 

= \atj(x,t)\ 

x [c 3 (x,t)e-^°^ + c](x,t)e^°^](19) 

Here we have the well known result that by properly ad- 
justing the local oscillator phase <j>j(x,t) one can mea- 
sure a particular quadrature component of the field. 

We now consider the correlation between two particu- 
lar quadrature components of the signal and idler fields, 
Z\ and Z 2l integrated over a finite detection time Td 

p(x i,x 2 ) = J Td dti J Td dt2(Z 1 (x 1 ,t 1 )Z 2 (x2,t 2 )}(20a) 
= J T dt\ J Td dt 2 al(xi,ti)a 2 (x 2 ,t 2 ) 

x{ci(xi,ti)c 2 (x 2 ,t 2 )) +c.c. (20b) 

The last line follows from adopting the parametric ap- 
proximation as well as from the fact that the signal and 
idler are in the vacuum state at the crystal entrance, 
since under these assumptions (c\(xi,ti)c 2 (x 2 ,t 2 )) = 0. 
In most cases Td 3> T co h so the detector is too slow to 
follow the temporal dynamics of the system. This means 
that the measured intensities are in fact proportional to 
the time-integrated field quadratures selected by the LOs. 

In the SPWPA we may evaluate the correlation l|20b[) 
using the formalism outlined in Eqs. l|5|l-l|10[l. We will do 
this in the degenerate case where signal and idler have 
identical wavelengths Ai = A2 = A and hence also vac- 
uum wave numbers fci iV = &2,v = k v — 2t:/X. In order to 
evaluate the correlation (|20b|l we note that Eq. (|15fl pro- 
vides the Fresnel transformation from the exit of the crys- 
tal to the detection plane. The correlations at the crystal 
exit are the near field signal-idler correlation which can 
be found from (|14bl) as 

(a^^xuh^&M) = 
J g e _ j0(tl _ ta) J _^_ e i^- S2 ) G% ^ (21) 



where we have introduced the gain function 

G(q,Sl) = \G{q,Sl)\e^^ = U x {q, Q)V 2 (-q, -SI). 

(22) 

Appendix 1X1 shows an approximate quadratic expansion 
of the gain phase that will become useful in the following. 

We now use (|15fl and (|21|l to obtain a more useful no- 
tation for the signal-idler correlation at the detectors 

(c 1 (if lj i 1 ) Ca (f2,t 2 ))= J ^e-^-t?) 

x J dqhi(xi,-q)h 2 (x 2 ,q)G(q,Q), (23a) 

r dx 1 

hj(x, q) = J ^e-**- x hj(x, Sf ). (23b) 

By using Eq. I|23a() we may evaluate the signal-idler corre- 
lation (|20b|l for Td ^> T co h- The details of this calculation 
are put in Appendix El an d basically there are two lim- 
its to consider. One is the case of a pulsed LO where 
the LO duration time is much smaller than the coher- 
ence time, and the result is given by Eq. l|B3Jl . The other 
case is that of a continuous- wave (cw) LO, where the LO 
duration time is much larger than the coherence time, 
and the result is given by Eq. JB5]). As mentioned in the 
appendix we did not find substantial differences between 
the two cases, and we therefore decided to use the latter 
in the rest of this paper for the analytical calculations. 
For completeness we repeat Eq. (|B5ll . stating that for a 
cw LO the following expression holds 

Td 

p(xi,x 2 ) = — a\{x x )al(x 2 ) 
Zir 

x J dqhi(x 1 ,~q)h 2 (x 2 ,q)G(q,0) + c.c, (24) 

where otj(x) = a.j(x, t = 0). 

In the following two sections we treat the two choices 
for the measurement technique in the signal arm, namely 
point-like detectors and bucket detectors. In each section 
both analytical and numerical results are presented. 



IV. POINT-LIKE DETECTOR SETUP IN TEST 
ARM 

When point-like detectors are used in the test arm the 
setup of Fig. \Ql is used. The signal arm we keep fixed 
in the f-f setup with the object at the crystal exit. The 
idler arm is then set in either the f-f setup or the tele- 
scope setup. The information is extracted by fixing the 
signal detector position x\ and scanning the idler detec- 
tor position x 2 . This gives a simple, fixed setup of the 
test arm. 
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A. Analytical results 

The analytical results are based on the SPWPA, so 
evaluating il23b|) for the signal we use (|16af> to obtain 

htixuq) = (iA/)- 1 f obj (f 1 fc v // + ^ ) (25) 

where T bj (q) is the Fourier transform of T bj (x) . 

1. Retrieval of object far field 

The object far field can be retrieved from the correla- 
tions in the point-like detector case by properly adjust- 
ing the optical setup in the reference arm based on the 
f-f scheme. We therefore evaluate l|23b|) in the f-f case 
where by using l|17a|l we obtain 

h2,f(x 2 ,q) = 2Tr(i\f)- 1 6(q + x 2 k v /f). (26) 

By inserting l|25|) and (|26[1 in (|24|) the signal-idler quadra- 
ture correlation is 

2T d 

p f (xi,x 2 ) = -jjj^\G(-x 2 k v / f,0)ai{x 1 )a 2 {x2)\ 

x Re (f obj [(£i + x 2 )K/f]e 1 *'^ A ^2,7a) 
<&f(xi,x 2 ) = 4>G(-x 2 k v / f, 0) 

-<fe°(x 2 )-<f>\°{x 1 )+-K. (27b) 

This is the main result for the retrieval of the object far 
field. We now discuss how to optimize the system for 
obtaining the far-field image information. 

First of all we notice that we need to keep the spa- 
tial dependence of the LO moduli as constant as pos- 
sible, effectively setting a limit on their spatial waists. 
However, since the effect of the LO moduli is a sim- 
ple multiplication on the correlation it is easy to do a 
post-measurement correction of the correlations. This 
means that matching the spatial overlaps of the LO mod- 
uli and the PDC fields should be straightforward, while 
the phases of the LOs are more crucial, as discussed be- 
low. The same comments turn out to hold for the other 
setups discussed in the paper. 

Second of all, when x 2 is scanned the gain modu- 
lus comes into play. We choose a non-collinear phase- 
matching condition (see Sec. IIVB|) in which \G(q = 

X 2 K X //, 0)| is a plateau shaped function centered on 
q = —qc, with [2(j 

qc = \p2l c ql = p 2 (l/fci + l/k 2 )-\ (28) 

Additionally, it is roughly constant over a finite re- 
gion determined by the spatial bandwidth go defined by 
Eq. I|12a|) . This sets a limit to the imaging resolution 
since the higher order spatial frequency components are 
cut off by the gain. In Sec. IIV B 31 we will show how to 
circumvent this limitation. 



Finally, the most crucial point is to control the LO 
phases to observe the desired quadrature components of 
the object far field. It is seen that the real part of the 
object far field is obtained if we can make $/(a?i, x 2 ) = 0, 
while if we can make 4> f{xi,x 2 ) = —tt/2 the imaginary 
part of the object far field is obtained. This can be 
achieved by engineering 02 /(a^) to cancel the gain phase 
dependence on x 2 that appears in Eq. (|27b(l . As a good 
approximation we may take </>2y (2:2) constant because 
the gain phase is a slow function over a region determined 
by q : an d we may therefore evaluate it anywhere inside 
the gain region, e.g. at the gain center —x 2 k v /f = —qc- 
Thus, to see the quadrature corresponding to the real 
part we should choose 

4>2°f = fai-fa 0) - 4>\°{x x ) + ir. (29) 

The quadrature component corresponding to the imagi- 
nary part of the object far field object may consequently 
be observed by adding tt/2 to this value. 

We can improve this result by using that, as shown 
in App. El the gain phase can be approximated with a 
quadratic expansion in q, i.e. 4>c(q,0) ~ + $~Q q ■ 

q + ^g^|g| 2 [see Eq. <|A3a|) ]. The linear term can 
be compensated by making the idler LO a tilted wave 
instead of a plane wave, i.e. 

$° f (x 2 ) = + ■ x 2 , (30) 

where ip 2 9 is a controllable reference phase and the wave 
number q^fj may be applied to the LO by using a grating. 
We should therefore choose the following value of the LO 
wave number 

92° = = P2lc* g k v /fe x . (31) 

The quadratic term 4> G q may be cancelled by shift- 
ing the focusing plane of the idler f-f system Az away 
from the crystal plane P [20|. since this gives (in 
Fourier space) a Fresnel transformation contribution of 
exp(— iAz|gj 2 /2fc v ) on the gain function in Eq. 121fl . 
Thus, using such a shift, the phase is approximately 

cPgaM 0) - 4 0) + 4 1} ■ Q + (4 2) - A^/2fc v )kf , 

(2) 

and by setting Az/2fc v = <Pq the quadratic phase term 
is exactly cancelled. So we must image a plane inside the 
crystal since this choice from Eq. (|A3d|) implies 

Az = -(1/m + l/n 2 )tf s Z c , (32) 

with *5> g defined in Eq. (|A2|) . Now all we need to do 
is to find the overall reference phase. To observe the 
quadrature component corresponding to the real part of 
the object far field we should choose 

^^("'-^(fO+i. (33) 

with (j^Q given by Eq. l|A3b|l . This result, the tilted wave 
LO (|31|) and the shift of the focusing plane (|32|l represent 
the optimized scheme for accessing the far-field image. 
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2. Retrieval of object near field 

The object near field can be observed by using the 
telescope setup in the reference arm. The idler kernel is 
then found from (|17c|) and (j23b(l as 

h 2 , T (x 2 ,q) = (2 7 r)- 1 e-^ 2 . (34) 
Using this in (|24|l we get 

x J ^e-V-^fouffiGfaky/f - q, 0) + c.c. (35) 

The correlation is therefore an integral over the gain func- 
tion and the object far field. Since the gain modulus is a 
roughly plateau-shaped function in a region determined 
by go that acts as a cut-off of the higher values of q, and 
the gain phase dependence is slow within this region, the 
gain can as a good approximation be pulled out of the 
integration as a constant evaluated at the center of max- 
imum gain q = —qc- The remaining integral is merely 
the inverse Fourier transform of the object far field, i.e. 
simply the object near field and we have 

,.. - x T d \G{-qc,0)ai(xi)a 2 (x 2 )\ 
p T {xi,x 2 ) ^ — 

7tA/ 

xRc[T obj (x 2 )e 4 *^ lA >], (36a) 

§ T {xi,X 2 ) EE -X 2 ■ Xikv/ f + (j>G(-qb,0) 

-#°(2i)-$£(* 2 )-7r/2. (36b) 

This is the main result for the retrieval of the object near 
field. We should stress that it is only an approximate 
result that neglected a cut-off of the high-frequency com- 
ponents of the object Fourier transform. Thus, image 
information that relies on large spatial frequency compo- 
nents will not be reproduced by the correlations. 

If we can make <^t{x\,x 2 ) constant and equal zero or 
— 7r/2, the real or imaginary part of the near-field object 
is observed. The linear contribution in x 2 can be can- 
celled by making the idler LO a tilted wave in the spirit 
of Eq. PU)1. i.e. 4> 2 ^(x 2 ) = ip^ + gf §. • x 2 and choosing 

<k% = -XiK/f. (37) 

In order to observe the real part of the object near field 
the overall constant reference phase must be chosen to 

V> 2 L ° = M-gb,Q) ~ #°(2t) ~ t/2. (38) 

We may optimize the result (|36a|l by compensating the 
quadratic term of the gain phase by taking the imaging 
plane inside the crystal by the amount (|32|l . in the same 
way as discussed in the previous section. In this case the 
idler reference phase should instead be 

^2° = <4 0) - + *i • 4 1)fc v// - it/2. (39) 



with <j>Q and <$q given in Eq. (|A3|) . 

This concludes the analytical part when point-like de- 
tectors are used for measuring the signal quadrature com- 
ponents. We emphasize that the retrieval of information 
about both the object near field and object far field can 
be done by solely playing with the optical setup of the 
reference arm. In the following subsection we show spe- 
cific examples from numerical simulations that support 
these results. 



B. Numerical results 

1. General introduction 

The numerical simulations presented in this paper take 
into account the Gaussian profile as well as the pulsed 
character of the pump. They consisted of solving the 
equations J5j using the Wigner representation to write 
the stochastic equations (see [20| for more details). At 
each pump shot a stochastic field was generated at the 
crystal input with Gaussian statistics corresponding to 
the vacuum state. Quantum expectation values were ob- 
tained by averaging over pump shots. This simulates the 
quantum behavior of symmetrically ordered operators, 
and thus the correlations calculated are symmetrically 
ordered. The system we have chosen to investigate is 
the setup of a current experiment performed in Como 
[23, H3- Thus, we consider a BBO crystal cut for type II 
phase matching. The pump pulse is at Ao = 352 nm and 
has To = 1.5 ps. The crystal length is l c = 4 mm and the 
dimensionless gain parameter a p l c = 4.0. Together with 
a pump waist of wq — 660 /im it corresponds roughly to 
a pump pulse energy of 200 fi3. The crystal is collinearly 
phase matched at an angle 9 — 49.05° of the pump prop- 
agation direction with respect to the optical z axis |29| . 
We used 9 = 48.2° which gives Ao < and a gain that 
is plateau shaped over a large region centered on —qc as 
given by Eq. ■ 

The correlations typically converged after a couple of 
thousands repeated pump shots when no spatial infor- 
mation was recorded in the test beam (i.e. using cither 
bucket or point-like detectors there) . For this reason it is 
very time-consuming to do a full 3+1D simulation (xyt 
propagated along z). Instead we chose to model a 2+1D 
system: a single transverse dimension x along the walk- 
off direction and the temporal dimension t with propa- 
gation along the crystal direction z. We will also show 
simulations where in addition to averaging over shots, a 
spatial average over x\ is also performed. In the 2+ ID 
simulations this reduces the number of shots needed by 
an order of magnitude, while in 3+ ID it is possible to 
average over both transverse directions and a diffraction 
image can be observed even with only a few shots. Thus 
this method allows us to calculate the correlations even 
from a full 3+ ID simulation. 

The integration along the crystal direction was done 
in N z = 200 steps. In the 2+1D simulations the art-grid 
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was chosen to N x — 512 and N t = 32, while for the 3+1D 
simulations we chose N x = N y = 256 and N t = 32. (x, t) 
were scaled to x co h and r co h, respectively, as given by 
Eq. 113|) . while z was scaled to the crystal length l c . 

We chose to investigate the quantitative behavior of 
the imaging system with a transmitting double slit as 
the object, defined in one transverse dimension as 

T (x) _ { 0, \ x -8 a \>zp, \x-5 x \<^ 

where d is the distance between the slit centers, a is the 
slit width, and S x is the object shift from origin. This ob- 
ject only introduces amplitude modulations and does not 
alter phase information. We take S x ^ in the numerics 
to compensate for walk-off effects, so we must choose it 
properly in order for each slit to see the same gain. The 
analytical Fourier transform is 

f bj(g) = ay27^e-^«cos(qd/2)sinc(qa/2), (41) 

where sinc(a;) = sin(x)/a;. The shift S x alone introduces 
a nonzero imaginary part. Additionally we will present 
results using a pure phase double slit [T(| 

TM _{-1, \x-5 x \>^,\x-S x \<^ 

This object does not alter amplitude information. 

In homodyne detection the LO is typically taken from 
the same source that created the pump. One would first 
generate the pump through second-harmonic generation 
(SHG), and before the laser light at the fundamental fre- 
quency enters the SHG crystal some of the power is taken 
out using a beam splitter, and this field is then used for 
the LOs. The second harmonic at the exit of the SHG 
process is then the pump for the PDC process. There- 
fore the LOs have more or less the same energy, shape 
and duration as the pump pulse entering the PDC setup. 
Through numerical simulations with Gaussian LOs we 
confirmed the results indicated in the analytical sections, 
namely that the effect of a Gaussian shaped LO appears 
trivially as a multiplication [see Eq. I|27a|l and l|36a|) ]. 
Since broader LOs can easily be achieved experimentally 
using a beam expander, we decided to keep the LOs as 
plane waves in space. Also this allows to base the inter- 
pretation of the results on the physics behind the PDC 
and keeping in mind that the LOs will only change the 
result with a multiplication of a Gaussian function. The 
temporal duration was chosen identical to the pump (1. 5 
ps), which ensures a good overlap between the pulses [26| . 

2. Averaging over shots 

In this section we present numerical results from aver- 
aging over repeated shots of the pump pulse, and where 
the optimization steps discussed in the analytical section 
were carried out (i.e. cancelling the quadratic and linear 
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FIG. 2: The reconstruction of the object far-field distribution 
using the f-f setup in the idler arm. Fhe numerical signal-idler 
correlations (open diamonds) were calculated from averaging 
over 2,000 repeated pump shots. In (a) and (b) ^2 ° was cho- 
sen as to observe the real and imaginary parts, respectively. 
The thin line shows T bj from Eq. I|41|l . while the thick line 
is Eq. <!27al as calculated in Mathematica. Slit parameters: 
a = 9 pixels, d — 33 pixels and 8 X = 23 pixels, corresponding 
to 34 fim, 123 fim and 86 /im, respectively, for / = 5 cm. 



phase terms of the gain). We also present semi- analytical 
calculations performed in Mathematica, which take the 
analytical formulas from Sec. IIV Al and carry out the op- 
timization steps and in addition take explicitly into ac- 
count the finite gain. 

The fixed position of the signal detectors X\ must 
be chosen with care. In the far-field case we must en- 
sure that in Eq. (|27a|l the object diffraction pattern 
TobjK^i + X2)k v /f) is centered where the gain modu- 
lus \G(— X2k v /f, 0)| has its maximum (which it has at 
—Xihv/f = —qc)- Thus we must choose Xxk v / f — 
—qc — —8.69. As for the near-field setup Eq. (|35() dic- 
tates that also here this position should be used to center 
the gain over the origin of Tobj. Thus, the test arm de- 
tector position is unchanged as we pass from the far field 
to the near field. 

Let us first take the case of an f-f setup in the idler 
arm. Figure [21 displays the two quadrature components 
corresponding to the real and imaginary parts, and this 
confirms that with this setup we can reconstruct the full 
phase and amplitude information of the object diffrac- 
tion pattern from the signal-idler correlations. The open 
diamonds in the figure are results from the numerical 
simulation. The results are compared to the analytical 
result (|41|) shown with a thin line, as well as a Mathemat- 
ica calculation (thick line) that effectively plots Eq. (|27all . 
Compared to the analytical object distribution, the nu- 
merical simulation shows that only the central part of 
the image is reconstructed since the high frequency com- 
ponents die out. This is a consequence of the finite gain 
bandwidth, as it was predicted by Eq. I|27afl : the Math- 
ematica calculation is based on this equation and it al- 
most coincides with the numerical results. An additional 
conclusion can be made from this agreement: since the 
Mathematica result assumes that the pump is a station- 
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(b) Imaginary part 
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FIG. 3: The reconstruction of the object near-field distribu- 
tion using the telescope setup in the idler arm. The notation 
and parameters are the same as in Fig. |5] except the thin line 
here shows T bj from Eq. 1401 while the thick line is Eq. 135H 
as calculated in Mathematica. x co h = l/f/o- 



ary plane wave this implies that the Gaussian profile in 
space and time of the pump pulse in the numerics does 
not affect the correlations significantly. This is because 
the object is in this example very localized so the signal 
field impinging on it is almost constant. At the end of the 
section we shall see an example, the pure phase object, 
where this is not the case. 

In Fig. |3 we show that by merely exchanging the f- 
f setup in the reference arm to the telescope setup, we 
can reconstruct the object near-field distribution from 
the signal-idler correlations. Fig. [31(a) shows the real part 
which follows the shape of the analytical double slit (|4l)H , 
while the almost zero imaginary part shown in (b) con- 
firms that the object does not alter phase information. 
The analytical SPWPA result using Mathematica is cal- 
culated from Eq. I|35[) [note that it is the integral over q 
that is calculated analytically and not the approximated 
version given by Eq. il36a[) ] . In comparison with the an- 
alytical object (thin line) only the basic shape can be 
recognized from the numerics, while as before the nu- 
merics agree very well with the analytical Mathematica 
result. This again indicates that the finite spatial band- 
width of the gain plays a decisive role in reconstructing 
the object. Actually, the results from the near-field and 
the far-field scenarios are linked: the far-field showed a 
cut-off at high frequency components due to finite gain. 
Exactly this cut-off in Fourier space is what gives the 
smooth character of the near-field image: as predicted 
by Eq. (|35|l the correlation is an integral over q of the 
gain multiplied with the object far-field, and because the 
finite gain bandwidth acts as a low-pass filter then it is 
not possible to reconstruct the sharp edges of the double 
slit as this would require high-frequency spatial Fourier 
components. 

The ghost-imaging schemes distinguish themselves 
from the Hanbury-Brown-Twiss (HBT) type schemes 
|30| in being able to pertain phase information of the im- 
age. In the HBT technique the object is placed in both 
of the correlated beams. When operated with thermal 
light the HBT technique can give information about the 
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FIG. 4: The reconstruction of a pure phase double slit in (a)- 
(b) the f-f setup case, and (c)-(d) the telescope case. The no- 
tation and parameters are the same as in Figs.|5|and|S| except 
10 4 shots are used. Note that the Mathematica results (thick) 
have been adjusted to take into account the finite shape of the 
signal field impinging on the object. 



Fourier transform of |T bj| 2 , while when operated with 
PDC beams it can give information about the Fourier 
transform of T^ b j [9l.ll2l|. Thus, the phase information is 
either lost or altered in the HBT scheme. Instead, in the 
ghost-imaging schemes the object is only placed in one 
of the beams, and from this one may instead reconstruct 
the Fourier transform of T bj . This holds both for PDC 
beams as well as for thermal (or thermal-like) beams ■ 

As a specific demonstration we used the pure phase 
double slit ^ as an object, since for this object the HBT 
scheme will not be able to provide information about the 
diffraction pattern. Figs. 0] (a) and (b) show the far-field 
reconstruction with the f-f system in the idler arm, and 
both the real and imaginary parts follow closely the ana- 
lytical Fourier transform (thin line). Thus, the scheme is 
able to reconstruct the Fourier transform of a pure phase 
object. However, there is a disagreement for the central 
peak in (a), corresponding to the average value of the 
real part of the near field: the analytical value actually 
goes out of the plot range shown. This is due to the fact 
that ideally the Fourier transform of the chosen phase 
object has a <5-like behavior here because all the photons 
are transmitted, whereas in the numerics the signal field 
impinging on the object has a finite extension. We ac- 
commodated for this in the Mathematica calculation by 
multiplying the phase double slit with a Gaussian fitted 
to the near-field profile of the signal in the numerics, and 
as can be seen (thick line) the agreement is now very 
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good. In Figs.|U(c) and (d) the near-field reconstruction 
is shown with the telescope setup in the idler arm, and 
as evidenced in the real part (c) the object image can 
be successfully reconstructed (the imaginary part is not 
shown since it was roughly zero). Also here the Math- 
cmatica calculations include the same finite size of the 
signal field and match very well with the numerics. We 
should stress that the "blurred" shape of the slits is due 
to the finite gain bandwidth (as in Fig.|3J, while the devi- 
ation away from the slits from the analytical result [thin 
line, Eq. (|42|) ] is due to the Gaussian profile of the signal 
[exactly this latter behavior gives rise to a lower central 
peak in (a)]. Finally, we also show the modulus of the re- 
constructed near field (d), and we do this to stress that it 
actually reveals some features of the double slit. This is 
a consequence of the finite gain bandwidth of the source 
which turns the signal-idler fields from being completely 
incoherent (this is the ideal result in the case of infinite 
bandwidth) to being partially coherent. Note that in the 
numerics we had to use more shots than before, otherwise 
the behavior away from the slits in the near-field case was 
not too regular. This is a consequence of the more noisy 
statistics there, as these parts have much lower intensi- 
ties. 



3. The spatial average technique 

In this section we describe how to improve the speed of 
the correlation convergence as well as to obtain a much 
larger imaging bandwidth. The idea is to average the cor- 
relations not only over repeated shots but also over the 
position of the signal detector in a certain way. This im- 
plies either using a scanning point-like detector setup or 
an array of detectors in the test arm. This abandons the 
idea about keeping the test arm setup as simple as pos- 
sible, but the benefits of having an increased bandwidth 
and convergence rate should not be neglected. The tech- 
nique works for reconstructing the far field, while for the 
near-field reconstruction this technique cannot be applied 
for several reasons. Most importantly, the final form of 
the idler correlation i|36al) when averaged over x\ gives 
the same result as using a bucket detector, and we al- 
ready argued in Sec. IIIII through the unfolded Klyshko 
picture that such a setup does not give any information 
of the object. Another factor is that $t(xi, x 2 ) depends 
linearly onfi, see Eq. (J36b(l and (|3"7Jl . Thus, it is impos- 
sible to engineer the phases to observe a given quadrature 
as fi is varied. 

For the f-f case, however, the technique is very powerful 
and it works as follows. Let us take the case where we 
want to see the real part, and assume we have performed 
all the optimization steps described in Sec. 11V A II The 
correlation is then 

Pf(xi,x 2 ) oc \G(~x 2 k v / f,0)ai{xi)a 2 (x 2 )\ 

xRe(f obj [(£i +z 2 )fcv//]) ■ (43) 



A change of coordinate system x = x\ + x 2 and averaging 
over x\ a,s pf^ 1 (x) = J dxiPf(xx,x) gives 



Pf,3 x (x*) oc Re f obj (xk v /f) 
x / dxi\G[(—x + xi)k v /f, 0}ai(xi)a 2 (x - x\)\ 



= Rc 



T bj{xk v /f) 



x J dt\G{Zk v /f,0} ai (Z + x)a 2 (-O\, (44) 

where in 144|) we have made a change of integration vari- 
able £ = x\ — x. Assuming that |ai(£ + x)\ is slowly 
varying with respect to £ in the gain region for any x 
inside the object diffraction pattern, we obtain 



p fySl (x) (xRe f ohi (xk v /f) \ai(x)\ 
x / dt\G& v /f,0}a 2 (-t)\, 



(45) 



The final integral is just a constant and does not depend 
on x, and hence as x is scanned the object far field is ob- 
served multiplied by the modulus of the signal LO. Thus, 
whereas with x\ fixed the imaging system had a finite 
bandwidth, now the bandwidth is effectively only limited 
by the shape of the signal LO (which consequently must 
be taken as broader than the far-field diffraction pat- 
tern). Therefore, if we assume that |ai(a?)| has a wide 
enough profile there is no cut-off of the spatial Fourier 
frequency components. The possibility of getting rid of 
the gain cut-off is is relying on the object only being lo- 
cated in one arm so the gain is a function of x 2 only 
while the object far field is a function of x\ + x 2 . In 
contrast, if the object is located in both arms, as in the 
HBT scheme, this would not be possible. Notice that 
this procedure in practice does not amount to merely in- 
tegrating over x\ (as when a bucket detector is used). 
Rather, one should move to position x\ of the signal de- 
tector while moving together the position x 2 of the idler 
detector as to keep x — x\ + x 2 constant. This corre- 
sponds to a spatial convolution between the signal and 
idler quadratures. To see that, consider the correlation of 
the measured quadratures Zj, and from Eq. (|20a|l using 
the substitution x = x\ + x 2 we have 



dt\ I dt 2 



x / dx 1 {Z 1 {x 1 ,t 1 )Z 2 (x' - x 1 ,t 2 )) (46) 



In the numerical simulations this convolution is rapidly 
calculated using the fast Fourier transform method. 

We remark that the bandwidth can also be improved 
without applying the spatial average. Namely, if instead 
x 2 is kept fixed and x\ is scanned, according to Eq. I|27a|) 
the gain cut-off is not present, while the diffraction pat- 
tern still emerges from the correlations. We performed 
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FIG. 5: The same as shown in Fig. except the average is 
performed over 200 pump shots as well as over xi. (a) and (b) 
show the real and imaginary parts of the reconstructed object 
far-field. Notice that the full numerical grid is shown, (c) and 
(d) show the real and imaginary parts of the inverse Fourier 
transform of the reconstructed far-field correlation function. 
The thin lines display the analytical results of Eq. (I4H for the 
far field and 14011 for the near field. 



numerical simulations that showed that this technique 
in fact gives an unlimited bandwidth as when a spatial 
average is performed over x\. However, since no spa- 
tial average is performed the convergence rate is as slow 
as when X\ is kept fixed. If an average is done over x-i 
we obtain the same result as the spatial average over x\. 
For a more complete discussion on extending the imaging 
bandwidth see Ref. |3l|. 

With the same system setup as in Fig. [21 the result 
for averaging over all x\ is shown in Fig. [S] First of all, 
we note the excellent agreement between the numerics 
and the analytical far-field pattern, and also that there 
is no frequency cut-off: the imaging bandwidth has be- 
come practically infinite. Moreover, only 200 shots were 
needed to obtain the same degree of convergence as in 
Fig. [21 which is an order of magnitude faster. Another 
interesting point about homodyne detection is that by 
measuring both quadratures in the the far-field distribu- 
tion, we may reconstruct the complete near-field object 
distribution from this information by using the inverse 
Fourier transform. We have done this for the data in 
Fig. (a) and (b) and the result is shown in (c) and (d). 
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FIG. 6: The same setup as used in Fig. 0] except now the 
average is also done over x\. The far-field quadrature corre- 
lations are calculated from 1,000 repeated pump shots, and 
(a) and (b) show the real part and modulus, respectively, of 
the inverse Fourier transform of the correlation. 



The real part, (c), of the inverse Fourier transform fol- 
lows T bj very precisely. This is because we now have 
access to many more high frequency components com- 
pared to the case shown in Fig. Thus, as the far-field 
imaging bandwidth is increased, the near-field resolution 
improves. It is instructive to note that in the absence of 
the spatial average the inverse Fourier transform of the 
far-field correlations shown in Fig. would give exactly 
the result reported with the telescope setup in Fig. [31 
This again underlines the strong link between the near- 
field and the far-field measurements when all the phase 
information is intact. 

Returning to the phase object of Fig. 0] we saw that 
with the signal detector position fixed the finite gain 
bandwidth makes the imaging system partially coherent 
and some phase information is present in the modulus of 
the correlations. For the phase double slit we repeated 
the simulations using the spatial average technique to re- 
construct the far-field distribution. The inverse Fourier 
transform of the correlation is shown in Fig. and the 
real part now shows excellent agreement with the ana- 
lytical phase double slit (thin line). Most importantly, 
Fig. Et> showing the modulus now does not reveal any 
information about the double slit. The extended imag- 
ing bandwidth achieved with the spatial average makes 
the imaging system truly incoherent and hence no phase 
information is transferred to the modulus. 

The spatial average over x\ works even better in 
two transverse dimensions because there are many more 
points to average over. In Fig.0we show two cases where 
different objects were used: (a) a simple amplitude trans- 
mission mask with the letters "INFM", and (b) a more 
complicated amplitude transmission mask showing a pic- 
ture of a wolf. The modulus of the inverse Fourier trans- 
form of the far-field correlation is shown for different 
number of shots. Evidently the simple mask (a) con- 
verges faster than the more complicated mask (b), but 
nevertheless in both cases a good, sharp image is ob- 
tained after relatively few shot repetitions. After addi- 
tional averaging over shots the irregularities disappear, 
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FIG. 7: The modulus of the inverse Fourier transform of the 
far-field correlation function as obtained in an f-f setup with 
two different objects (a) and (b). The correlations are cal- 
culated from a full 3+ ID simulation, and averaging over the 
entire xi plane and over the number of repeated pump shots 
shown on the right. 



as indicated in the last shots. 

We now evaluate the qualitative difference of the near- 
field resolution between the fixed detector case and the 
spatial average technique. Fig. |5] (a) displays the mod- 
ulus of the near-field correlations obtained by using the 
telescope setup and a fixed signal detector position, and 
Fig. [S] (b) displays the modulus of the inverse Fourier 
transform of the far-field correlations obtained using the 
f-f setup and applying the spatial average technique. The 
blurred character in the fixed detector case that was also 
present in one transverse dimension (see e.g. Fig. [2J is 
now very obvious, and does not improve much even if 
more averages are performed since it is a consequence of 
the finite bandwidth as mentioned before. In contrast 
the much improved far-field bandwidth of the spatial av- 
erage technique increases the near-field resolution so the 
observed near-field image completely sharp. Note that 
the simulations in Fig. |H1 have been performed neglecting 
the time dimension, without which a dramatic extension 
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FIG. 8: Qualitative comparison of the resolution in (a) the 
fixed detector case and (b) the spatial average case, (a) uses 
the telescope setup in the idler arm to obtain the near-field 
correlation, and the modulus is shown, (b) uses the f-f setup 
in the idler arm to obtain the far-field correlation, and what 
is shown is the modulus of the inverse Fourier transform of 
the correlation. The simulations neglect the time dimension 
and 5,000 shots are used for the averages. 



in the amount of computational time would occur. Such 
an approximation corresponds to using a narrow interfer- 
ence filter, and the approach was justified by comparing 
the correlations at a fewer number of shots with simula- 
tions including also time. 

We have in this section shown that using point-like de- 
tectors in the test arm allows to reconstruct both ampli- 
tude and phase information about an image. The near- 
field and far-field object distributions could be obtained 
by only changing the optical setup of the reference arm. 
The spatial imaging bandwidth for a fixed test detector 
position is determined by the source bandwidth, and we 
showed that this in turn determines the near-field res- 
olution. Using a spatial average technique improves the 
imaging bandwidth dramatically in the far-field case, and 
also leads to a much faster information retrieval. 



V. BUCKET DETECTOR SETUP IN TEST 
ARM 

We now turn to the case where bucket detectors are 
used in the test arm. One motivation is to show that by 
using homodyne detection together with bucket detec- 
tors, phase-only objects can be imaged. If an intensity 
detection scheme is used this is not possible, as shown 
in Ref. 0] for the two-photon coincidence imaging case 
(although it is pos sible using a point-like detector in the 
object arm pi llO|V This is a general result that also 
holds in the high-gain regi me (it can easily be derived 
from the results of Ref.0, Another motivation is 

to see the imaging capabilities of this setup since it tech- 
nically seems simpler than the point-like detector setup. 
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A. Analytical results 

As argued in Sec. 11111 when bucket detectors are used 
in the test arm it is necessary to keep the object placed 
in the measurement plane. Besides that the lens setups 
in the arms are arbitrary, but we decided to keep the f-f 
setup in the signal arm, as well as the idler setup in cither 
f-f or the telescope configuration. The test-arm set-up is 
shown in Fig. The signal kernel (|16b|l is transformed 
according to i|23bjl as 

h\{x u q) = ^j5(q + x 1 k v /f)T ohi (x 1 ). (47) 

The idler kernels remain unchanged. The bucket detec- 
tors in the test arm effectively corresponds to an integra- 
tion over X\ of the signal-idler quadrature correlation as 
p^(x 2 ) = J dxip\(xi,x 2 ). 

1. Retrieval of object near field 

For an f-f setup in the reference arm the idler kernel is 
given by Eq. I|26|) and from Eq. (|24ll we therefore obtain 

p h f {xi,x 2 ) = — S(xi + X2)\a 1 (x 1 )a 2 (x 2 )G(xik v / f,0)\ 

x RelTobj^Oe 1 */^^)], (48) 
The measured signal-idler correlation is consequently 

rp 

p° f (x 2 ) = — \ai(-X2)a 2 (x2)G(-x 2 k v / f,0)\ 

x Re[T obj (-x 2 )e*/^)], (49a) 

$ b f (x 2 ) = <M-f 2 fc v //,o) 

-<ft°(-g 2 )-$°(2 2 )+n. (49b) 

This is the main result for retrieving the near field in the 
bucket detector case. 

The idler LO phase may then be engineered to observe 
the desired quadrature in much the same way as shown 
in in Sec. lIV ATI As a first approximation the gain phase 
dependence is neglected, so to observe the real part we 
must have ^(x 2 ) = 0, implying we should choose 

rij = M-fc,0) - #°(-* 2 ) + 7T. (50) 

We may optimize this result by setting the focus plane 
of the f-f setup in the idler arm inside the crystal with the 
amount given by Eq. I|32|l as to compensate the quadratic 
term of the gain phase. Additionally, by making the idler 
LO a tilted wave with the wave number given by Eq. I|31(l 
we compensate for the linear phase term. With these 
optimizations, the real part can be observed by choosing 

V2 L ,? = 4 0) -^ O (-^) + 7r. (51) 

We note that the results show a dependence of the 
signal LO as a\{—x 2 ). Therefore, as the idler pixels are 
scanned it is crucial that the signal LO does not vary 
substantially. Another point is the scale of the coordinate 
system, but we will return to this in the next subsection. 



2. Retrieval of object far field 

For a telescope setup in the reference arm the idler 
kernel is given by (|34J) and thus using Eq. (|24J) we have 

Py(xi,x 2 ) = 2ni\f a * 1 ) a 2 (x 2 )G(xik v / /, 0) 

xTobjOriK^-^^+c.c. (52) 
Thus, the measured correlation is 
T 

Pt{x 2 ) = -^a* 2 (x 2 ) 

X / ^■e- ai - 5f! ' fc -/'G(£ifc v //,0)aI(^i)T obj (£i) 

+c.c. (53) 

We assume now that the plateau-shaped gain and the 
modulus of the signal LO vary slowly with respect to the 
object, so they can be taken out of the integral, evaluated 
at the position of maximum gain x,\k w j f = —qc- This 
approximation yields 

p h T {x 2 ) ~ ^|G(-gb,0)ai(-gb//fcv)a a (52)| 

xRc[T obj (f 2 fc v // + q^°)e^ s % (54a) 
$£(x 2 ) = M-fc,0) - - 4>\° T {X2) + vr/2. (54b) 

This is the main result showing that the object far field 
can be reconstructed in the bucket detector case. 

We can observe the real of the far-field of the object if 
$^(x 2 ) = 0, which in the case of (|54b(l implies 

$° T (X2) = M-fc, 0) - + tt/2. (55) 

As in the other cases the result ()54all can be optimized 
by compensating for the quadratic gain phase term be- 
fore the integration by taking the imaging plane of the 
telescope setup inside the crystal with the amount given 
by Eq. I|32l) . In that case we should choose the following 
reference phase for the idler LO to observe the real part 

02,t(^) = 4 O) -^ O + 7t/2. (56) 

In connection with optimizing the imaging performance, 
it is apparent that if the signal LO is a tilted wave 
then the origin of the reconstructed diffraction pattern 
changes. This can be used to compensate for undesired 
walk-off effects by centering the reconstructed image in 
the place where the idler (which in this setup is in the 
near field) has its maximum. Additionally, making the 
idler LO a tilted wave can be used to compensate for 
unimportant oscillations in the quadrature components 
that comes from shifting the object near field from origin 
[see in that connection Eq. (|41|l ] . 

Let us come back to the scaling of the image infor- 
mation. The particular setup we chose (having the ob- 
ject placed in the far field created by an f-f lens system, 
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FIG. 9: The reconstruction of the object near-field distribu- 
tion of a pure phase object with bucket detectors in the test 
arm and using the f-f setup in the reference arm. The numeri- 
cal correlations (open diamonds) were calculated from averag- 
ing over 10 4 pump shots. The thin line is the analytical phase 
double slit, while the thick line is calculated in Mathematica 
on basis of Eq. I49al . Slit parameters in pixels as Fig. cor- 
responding for / = 5 cm to a — 166 fim, d — 610 /im and 
5x = —2.9 mm. 



just before the 
nal field impinj 

x f = fqo/K- 
mulas Eqs. 
ject was placed 
the scale x co h = 
and A = 704 
Xcoh = 17 /im. 
is possible with 
(i.e. larger) lenj 



bucket detectors) implies that the sig- 
ing on the object changes on the scale 
This is evident in the analytical for- 
and H53|) . In contrast, when the ob- 
in the near field the field changed on 
= 1 /go- As an example, for / = 5 cm 
nm, Xf — 338 /im while in contrast 
The consequence is that while imaging 
bucket detectors, it occurs on a different 
;th scale than with point-like detectors. 



B. Numerical results 

The numerical simulations in the bucket detector case 
were done exactly as discussed in Sec. IIVBI except now 
the object is placed after the f-f setup in the test arm. We 
decided to keep the object structure used in the point- 
like detector case, i.e. exactly the same number of pixels 
between the slits as well as for the slit aperture. As 
discussed above, this means that the slit dimensions are 
determined by the larger scale xj. Moreover, the pure 
phase double slit in one transverse dimension is used as 
an object in order to demonstrate the homodyne detec- 
tion scheme's capabilities. Finally, we mention that in 
the numerics all the optimization procedures discussed 
in the previous sections were carried out. This includes 
in the telescope case to use a tilted wave signal LO to cen- 
ter the reconstructed diffraction pattern over the peak of 
the idler near field, as well making the idler LO a tilted 
wave as to cancel the oscillations in the diffraction pat- 
tern originating from centering the object over the signal 
far-field plateau. These latter two optimizations are not 
crucial for the imaging performance. 

Fig. [5] shows the numerical simulation of the recon- 
struction of the object near field using the f-f setup in 
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FIG. 10: The same as shown in Fig. [5] except now using 
the telescope setup in the idler arm making it possible to 
reconstruct the object far-field distribution. The thick line is 
Eq. 153H calculated in Mathematica while taking into account 
the finite shape of the idler near-field distribution. 



the idler arm. The reconstructed real part in (a) follows 
the analytical double slit, but notice that the Mathemat- 
ica result (thick line) has a better resolution than the 
numerical results. This indicates that the finite shape of 
the pump worsens the resolution, and we will come back 
to this later. The imaginary part (b) is containing only 
little information about the slits. Away from the slits the 
correlations decay to zero because of the finite shape of 
the far-field gain. 

Fig. 1101 shows a numerical simulation similar to Fig. El 
but with the telescope setup in the idler arm allowing 
for the reconstruction of the object far-field distribution. 
The reconstructed real part in (a) follows within a cer- 
tain bandwidth the analytical result (thin line) . The cut- 
off is actually caused by the Gaussian shape of the idler 
near-field profile, which has been taken into account in 
the Mathematica calculation of the integral i|53|) ; without 
such a correction the Mathematica result would actually 
very closely follow the analytical result (thin line) with- 
out being cut off. This can be understood by inspecting 
Eq. 1 ■")/)! more closely, since we see that the action of the 
gain inside the integral is to provide a limit for the object 
extension. However, in contrast to the point-like detector 
case of Eq. I|27a|l where the gain eventually determines 
the far-field imaging bandwidth, in the bucket case it is 
not so: Once the object is located inside the gain, the 
reproduced diffraction pattern is very good. This is also 
in accord with what we saw in the semi-analytical repro- 
duction of the near field in Fig. [51 where the double slits 
appear very sharp indicating that the Fourier frequency 
bandwidth is large. A consequence of this is that if we 
tailor the pump to have a plane-wave shape, we should 
see practically no cut-off in the far field as well as a very 
sharp reproduction of the double slit in the near field. We 
checked this to be the case. Thus, also the bucket detec- 
tor case shows an inherent link between far-field band- 
width and near-field resolution. Finally, note that the 
imaginary part (b) is zero because we chose to eliminate 
the oscillations arising because of the offset of the object 
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from origin, making the Fourier transform purely real. 
Additionally, the central peak of the analytical function 
in (a) goes out of the range shown for the same reasons 
as in Fig. ^ 

In conclusion, a homodyne detection scheme with 
bucket detectors in the test arm is capable of observing 
both amplitude and phase distribution of an object. Both 
the near-field and the far-field distributions are accessible 
by only changing the optical setup in the reference arm 
and even phase-only objects can be imaged. 



VI. CONCLUSION 

We have shown how a homodyne detection technique 
can be used to get access to both modulus and phase in- 
formation in ghost-imaging schemes based on parametric 
down-conversion in the high-gain regime. It is necessary 
to measure both the signal and idler beams with individ- 
ual local oscillators, and we have shown analytically how 
to engineer the phases of these in order to retrieve both 
the object image (near field) and the object diffraction 
pattern (far field). The results were confirmed with nu- 
merical simulations taking into account a finite shape of 
the pump pulse. 

We showed that the resolution in the near field and the 
imaging bandwidth in the far field were inherently linked, 
and were restrained to the finite gain of the source. Thus 
in the far-field setup the gain cuts off the higher frequency 
components in the spatial Fourier domain. In the near- 
field setup the image is more blurred than the reference 
image, and we showed that this is exactly a consequence 
of the higher frequency components missing in the spatial 
Fourier domain because of the gain cut-off. 

The homodyne scheme allows us to have access to both 
phase and amplitude information, implying it is possible 
to pass from the far-field result to the near-field result 
by simply making an inverse Fourier transform. Thus 
by taking the far-field correlations containing the infor- 
mation about the diffraction pattern and performing an 
inverse Fourier transform gave the same result as measur- 
ing the near-field correlations containing the information 
about the object image, and vice versa. This also means 
that once we have measured e.g. the far field there is 
actually no need to also measure the near field. 

We showed that in the far-field setup a spatial average 
technique over the signal detector pixels leads to much 
faster convergence rate of the correlations (fewer pump- 
shot repetitions needed) as well as a hugely improved 
bandwidth of the diffraction pattern. An intuitive expla- 
nation of the method is that each far-field mode is un- 
correlated to the rest of the modes in the beam, but the 
strong mode correlation between the signal-idler beams 
means that for each signal mode an independent correla- 
tion can be measured containing the object information. 
The spatial average then corresponds to averaging over 
all the modes in the signal beam. This means that much 
fewer repeated pump shots of the laser are needed be- 



fore the correlations converge, and this truly exploits the 
possibility for parallel operations in spatially correlated 
beams. The increased bandwidth stands out as the true 
advantage of the technique, however, and it relies on the 
fact that as we change the test detector position a differ- 
ent part of diffraction pattern is amplified. Thus, as we 
average over all the positions we effectively also extract 
information about different parts of the diffraction pat- 
tern, resulting in an increased bandwidth. It is important 
to note that the technique also works when using direct 
intensity measurements of the fields instead of homodyne 
detection, both when using PDC beams (as will be dis- 
cussed elsewhere in a separate publication j3j|) and the 
thermal-like beams investigated in [l^. The technique 
works particularly well in the high-gain regime because 
the number of signal photons per mode is not small, but 
should also work in the low-gain regime. 

In the time domain the outcome of a homodyne de- 
tection is depending strongly on the temporal overlap 
between the local oscillator and the field to be measured. 
A question that had to be answered in this paper was 
therefore how important it is to ensure a proper overlap 
also in the spatial domain between the local oscillator 
and the field. The answer in the context of ghost imag- 
ing is that it is not crucial. In fact, the analytical results 
in the plane-wave pump approximation suggest that the 
near-field shape of the local oscillator is merely multiplied 
onto the final result. Even the numerical simulations tak- 
ing a Gaussian shape of the fields into account suggest 
the same. Another issue entirely is measurements where 
the quantum efficiency is important (in contrast to here), 
and in that case undoubtably the spatial overlap becomes 
more important. 

Using an intensity detection scheme the detection time 
must not be larger than the coherence time, otherwise 
the background term that contains no image informa- 
tion becomes dominating and the visibility of the image 
is dramatically decreased 0, 0, El- In contrast, the 
homodyne detection scheme is not restricted by this be- 
cause the correlation function is second order in the field 
operators implying there is no background term. Hence 
the detection time may be chosen much larger than the 
coherence time. Additionally, we have shown that the 
spatial average technique can give a much larger imaging 
bandwidth than the source itself. This suggests the pos- 
sibility of using an optical parametric oscillator (OPO) 
as a source, and using homodyne detection as a mea- 
surement protocol. The limited spatial bandwidth of 
the OPO with respect to PDC could then be circum- 
vented with the spatial average technique, and the homo- 
dyne measurements allow a cw operation of the imaging 
scheme: since there are no problems with image visibil- 
ity, the measurement time can be taken longer than the 
coherence time. Thus, one could effectively abandon the 
involved scheme of using a high-power laser pumping a 
PDC setup and the complicated task of overlapping LO 
fields with the signal-idler fields with temporal durations 
on a pico-second level. 
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APPENDIX A: THE APPROXIMATE FORM OF 
THE GAIN PHASE 



It is useful to obtain an approximate expression for the 
gain phase <pc from the expressions We note in this 
connection that Ai2(<f, ft) = A2i(— q, —ft) = A(q, 51) and 
therefore also T 12 (q, 51) = T 21 {-q, -ft) = T(q, 51). Thus, 
we obtain 



-AqI c + arctan 



A(q, ft) tanh[r(g, ft)l c 



2T(q, 51) 



-Aol c + arctan [ty g A(q, ft)l c 
-A l c + ^ g A(q,ft)l c . 



(Al) 



The approximation to Eq. ijAlf) is that the gain is large 
and that the modes inside the gain are phase matched 
(A(q, ft) ~ 0) making T(q, ft) ~ a p . We have also intro- 
duced a dimensionless phase parameter that is related to 
the dimensionless gain parameter o p l c 



tanh (a p l c 

2o~ p l c 



(A2) 



For a p l c —* (small gain) ty g 



1/2 while for a p l c 1 
(large gain) <3> ff — > 0. Now from Eq. i|10c|) it is evident 
that A(<f, ft) is quadratic in q and 51, so that the gain 
phase I|A1J) approximately can be written as 



<f> G (q, ft) ~ 4>f + $%] q • q + ^ q \q\ 



I„12 



4V + 4V 2 > ( A3a ) 



4 0) = A l c (-l + * g ), 
,(2) _ {ni+n 2 )lc^ g 

*G,q - 



2niri2k v ' 



\3X2 — (^1 k 2 )lc^g 



<& = [K+k'i)i^ g /2 



(A3b) 
(A3c) 

(A3d) 

(A3e) 
(A3f) 



where e x is a unit vector in the cc-direction. We shall use 
these expressions when the optimal phases of the LOs 
come into question. 



APPENDIX B: THE SIGNAL-IDLER 
CORRELATION 

Here we present the detailed calculations of the general 
expressions for the signal-idler correlation i|20b(l . 

By using Eq. I|23a|l we see that in the signal-idler cor- 
relation (|20b|) the evaluation of the temporal part only 
concerns the gain and the LOs, while the spatial im- 
age information comes out from the integral over q in 
Eq. (|23al) . To evaluate the temporal part of the correla- 
tion we rewrite (I20bl) as 



p(xi,x 2 ) = J dflB(xi,X2,ft) 
x Jdqhi(xi,— q)Ii2(x2,q)G(q, ft) + c.c, (Bla) 

B(^,x 2 ,ft) = T^JJ ^fal (x, ,fl')a* 2 (x 2 , 51") 

xsinc[(0 - fi')T d /2}smc[(ft + fl")T d /2], (Bib) 

Let us consider the case of a pulsed LO, i.e. evaluate 
Eq. i|Blb|) for the LO duration much smaller than r co h. 
Additionally taking into account that, as mentioned be- 
fore, it can be assumed that T d r co h we may use the 
approximation T d sinc[($! ± ft')T d /2] ~ 2TrS(fl ± ft'). In 
this limit we therefore obtain 



B{x u x 2 ,ft) = al(xi,ft)a 2 (x 2 ,fl), 
which means that the correlation is 



(B2) 



p(xi,x 2 ) = J dftal(xi,ft)a 2 (x 2 ,fi) 

x J dqhi(xi,-q)h 2 (x 2 ,q)G(q, ft) + c.c. (B3) 

In the case of a continuous wave LO, the LO duration 
is much longer that r co h- Thus, the LO corresponds to 
a quasi-monochromatic wave, i.e. ctj(0, ft) = 5(fl)aj(x). 
Using this form in Eq. (|Blb|l removes the integration on 
ft' and ft" . Then we use that for T d 3> T co h the remaining 
term T d [snic(ftT d /2)] 2 behaves like 2ir5(fl), implying 



B(x u x 2 ,ft) ~ 5{fi)a* 1 {x 1 )a* 2 (x 2 )T d /(2TT), 
which means that the correlation is 



(B4) 



p(xi,x 2 ) 



2n 



al(xi)a 2 (x 2 ) 



X / dqhi(xi, — q)h 2 (x 2 , q)G(q, 0) + c.c. 



(B5) 



It is clear that the pulsed case may have a larger 
gain than the monochromatic case because of the in- 
tegration over ft. However, a crucial point to achieve 
this goal lies in the minimization of the ^-dependence 
of the phase in the integrand of Eq. I|B3(1 . This will 
in turn maximize the integral regardless of the value 
of q. This can be done approximately by putting in 
a proper temporal delay between the LOs, effectively 
cancelling the first-order term i|A3e|) of the gain phase. 
The second-order term cannot easily be cancelled, but 
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we estimate that this term does not contribute much 
since we have a negligible contribution of this term of 
(k'{ + A£%n§/2 = (k'{ + fc£)/p2 - K) 2 ic] = Q (1(T 4 ) 
for the crystal setup we consider (see Sec. liVBTl for de- 
tails on the setup). Numerical simulations showed that 
a LO optimized in this way as expected enlarged the ef- 



fective spatial bandwidth. However, the improvement is 
not substantial and furthermore it requires temporally 
short pulses (on the order of 20% of the coherence time, 
or less) which are difficult to obtain experimentally. We 
therefore decided only to show analytical results in the 
limit of a continuous wave LO. 
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